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Steady-state fluorescence spectra of a two-level atom embedded in a three-dimensional photonic 
bandgap crystal and driven by a monochromatic classical electrical field is calculated theoretically 
for the first time as we know. The non-Markovian noises caused by the non-uniform distribution of 
photon density of states near the photonic bandgap are handled by a new approach in which the 
Liouville operator expansion is utilized to linearize the generalized optical Bloch equations. The 
fluorescence spectra are then directly solved by the linearized Bloch equations in the frequency 
domain. We find that if the atomic energy level is far from the bandgap, fluorescence spectra with 
MoUow's triplets are observed. However, when the atomic energy level is near the bandgap, the 
relative magnitude and the number of the fluorescence peaks are found to be varied according to 
the wavelength offset. 

PACS numbers: 32.50.-|-d, 42.70.Qs, 42.50.Pq 
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I. INTRODUCTION 

In recent years, with the advance of new fabrication 
technologies, it has become more feasible to actually uti- 
lize higher dimensional periodic dielectric structures (or 
especially the photonic bandgap crystals) P, for mod- 
ifying the properties of the photon states as well as the 
properties of the spontaneous emission. Like electrons 
in solid state crystals, photons are prohibited to prop- 
agate inside photonic bandgap crystals due to the lack 
of available photon states. In general, the photon den- 
sity of states (DOS) of photonic crystals has a highly 
non-uniform distribution near the bandgap, which is to- 
tally different from the case in free space. Such a non- 
uniform distribution of the DOS has been investigated by 
many authors and has provided a new and experimentally 
feasible platform for investigating photon-atom interac- 
tion. Many new phenomena such as photon-atom bound 
states 0, spectral splitting 0, quantum interference 
dark line effect 0|, phase control of spontaneous emis- 
sion 6] , transparency near band edge 'T'I , and single-atom 
switching Q have been theoretically discovered in the 
presence of the bandgap. From the Aulter-Townes spec- 
tra for atoms coupled to a photonic bandgap structure 
Q or equivalently a frequency-depended photon DOS 9J, 
the modification of the spontaneous emission caused by 
the environment [To| actually can be verified. However, 
all of the above studies only focused on the transient be- 
havior of the atom-photon interactions and to the best of 
our knowledge there is still no theoretical treatment on 
calculating the steady-state fluorescence spectra in pho- 
tonic bandgap crystals. 

To investigate this problem, the approach of the 
present paper is to treat the photon states of the photonic 
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crystal as the background reservoir and introduce non- 
Markovian noise operators caused by the non-uniform 
DOS distribution near the band edge. To model the non- 
uniform DOS distribution near the bandgap, an isotropic 
dispersion relation ^ and an anisotropic dispersion re- 
lation have been proposed to serve as a simple but 
qualitatively correct model for theoretical analyses. With 
the use of the three-dimensional anisotropic dispersion re- 
lation, we will investigate the steady-state properties of 
the resonance fluorescence spectra emitted by a two-level 
atom which is embedded in a photonic bandgap crystal 
and driven by a classical pumping light. Due to the non- 
Markovian nature of the atom-field interaction in this 
case, we can not directly utilize the Born-Mar kovian ap- 
proximation which is usually used in quantum optics for 
treating atom-photon interaction problems. To overcome 
this difficulty, we first derive the generalized optical Bloch 
equations without any approximation and then use the 
zero-order Liouville operator expansion to approximate 
the nonlinear terms. By solving the simplified linear 
equations in the Fourier domain directly we can calcu- 
late the stationary two-time correlation functions of the 
atomic operators as well as the spectral distribution of 
the resonance fluorescence. When the emission frequency 
of the atom is far from the bandedge, the triplet spectral 
shape is obtained, just as one will expect intuitively from 
the results first calculated by B. MoUow 12] for the free 
space case. But when we change the emission wavelength 
of the atom to be close to the band edge, both the num- 
ber of the peaks as well as their profiles are found to be 
varied depending on the wavelength offset between the 
atomic transition wavelength and the band edge. The 
details of these results will be reported in the rest of this 
paper. 

The paper is organized as follows: in section ^ we 
derive the generalized optical Bloch equations with noise 
operators caused by the surrounding reservoir. In section 
mil we use the anisotropic dispersion relation for model- 
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ing the photon DOS of the three-dimensional photonic 
bandgap structure and based on this model, we present 
the calculated fluorescence spectra in section llVl Finally, 
a brief conclusion is given in section Ivl and the first or- 
der Liouville operator expansion is carried out in the ap- 
pendix for accuracy checking. 



II. THEORETICAL MODEL 

With the help of rotation wave approximation, one can 
use the Jaynes-Cummings model to describe the atom- 
photon interaction ^|. For a single two-level atom em- 
bedded in a photonic bandgap crystal and driven by a 
monochromatic classical pump light, the photon states 
in the photonic bandgap crystal can be treated as the 
background reservoir. The Hamiltonian of the system 
can be written as: 



tors, we get the generalized Bloch equations as follows: 



k 

k 



(1) 



(3) 



(4) 



dt'Gc{t - t')a+{t')a^{t) + n+{t) 
-2 / dt'[G{t-t')a+{t)a-it') + Gcit-t')a+{t')a-{t)] 



Here A = lul ~ uja, and = oja — (jJk- The two func- 
tions, G{t) and Gc{t), are the memory functions due 
to atom-reservoir interaction and are defined as G{t) = 



a+e— *) Ek\9k\'e^^^^e{T), and G^r) ^ Ek\9k\'e 



e(r). 



Here 8(t) is the Heaviside step function. Moreover, 
the three noise operators n^{t), and are ex- 

pressed as follows: 



where the transition frequency of the atom and the fre- 
quency of the pumping light are denoted by uja and lol 
respectively, and ak are the creation and annihilation 
operators of the background field reservoir, fl is the Rabi- 
flopping frequency of the atom under the external pump- 
ing light and also represents the relative magnitude of the 
pumping light, a, = (|2)(2| - |1)(1|), a+ = |2)(1| = al 
are the usual Pauli matrices for a two-level atom, and 
gk is the atom- field coupling constant. Here we use the 
index k to label different photon states. The form of the 
coupling constant gk depends on the gauge one choose. 
Although the Lamb shifts calculated in non-relativistic 
quantum field theories are different with the p ■ A and 
f • E formula Q| , one can use the p ■ A formulation to get 
the correct form for the atom-field coupling if the 
electromagnetic field varies little over the spatial extent 
of the electronic wave function. Under the p ■ A formula- 
tion the coupling constant can be expressed as: 



5fe(d,ro) = gk = \d\uja 



1 



2h€oU!kV 



d-E*(fo) (2) 



Here we use notations |d| for the magnitude of the atomic 
dipole moment, d for the unit vector along the direction 
of the dipole moment, V for the volume of the quantiza- 
tion space, and eo for the Coulomb constant. 

The evolution equations for the atomic operators can 
be derived directly from the Hamiltonian in Eq. . By 
transferring the system to the rotating frame with the 
frequency uJa and eliminating the reservoir field opeva,- 



n^{t) =t"^gke'^>''<j,{t)ak{~x>) (6) 

k 

n+{t) = -^Y,gle~^^^'ai{-^)aAt) (7) 

k 

n,{t) = -2i^gfce*^'=V+(t)afc(-oo) (8) 

k 

+2z^5fee-*^^*at(-oo)a_(i) 



Supposing that the reservoir is in thermal equilibrium, 
then the mean and correlation of the reservoir field op- 
erators before interaction will be: 



{ak{-oo))F, = {a\{~(x,))R = 
{ak[~(X))ak'{-co))R = 
(4(-oo)4,(-oo))h = 
{a\{-oo)ak'{-'x>))R = UkSkk' 
{ak{-(X))a\, {-co)) R = [uk + l)5kk' 



(9) 
(10) 

(11) 
(12) 

(13) 



Here Uk is the mean quantum numbers of the reservoir 
modes under thermal equilibrium. 

Using the statistical characteristics of the reservoir 
field operators, it can be easily shown that the three 
noise operators ri_(i), ?T.+ (i), and nz{t) are zero mean 



3 



with their correlation functions given below: 

{n-{t))R = {n+{t))R^{n,{t))R = 
= {n+{t)n+{t'))n = 

k 
k 
k 

Since in general the correlation functions of these 
noise operators are not delta correlated at time (non- 
Markovian) , we cannot directly use the Born-Markovian 
approximation to solve the problem. One can see that the 
correlation functions depend not only on the photon den- 
sity of states, but also on the correlations of the atomic 
operators. Eqs. H3I5|) are called the generalized optical 
Bloch equations and will serve as the starting point for 
our further derivation. 
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FIG. 1: Amplitude and phase spectra of the memory function 
G{uj) with ujc = 100/3. The memory function is non- uniform 
around the bandedge and becomes pure imaginary inside the 
bandgap. 



III. MODELING OF THREE-DIMENSIONAL 
PHOTONIC CRYSTALS 

Due to the bandgap characteristics of photonic crys- 
tals, the spectral distribution of the photon DOS will 
have one or more discontinuities and is highly non- 
uniform near the bandedge. Although in general the 
DOS of photonic bandgap crystals will vary with the ge- 
ometrical structure and the dielectric constants of the 
materials they are made by, it is still possible to intro- 
duce some simple formula for approximately describing 
the DOS behavior near the bandgap. In the literature an 
isotropic dispersion relation Q and an anisotropic disper- 
sion relation jTll | have been proposed to serve as a simple 
but qualitatively correct model for theoretical analyses. 
However, according to the results from a full vectorial 
numerical calculation 16], the DOS near the bandgap 
for a three-dimensional photonic crystal increases from 
zero and behaves more like the anisotropic model. This 
is why we will adopt the anisotropic model in the fol- 
lowing derivation. From Eq. (0), one knows that the 
magnitude of the coupling constant is also dependent on 
the local electric field. Although in general one should 
apply this position-dependent coupling constant with the 
actual DOS to fully describe the photon-atom interaction 
within a photonic bandgap crystals, however for simplic- 
ity we will use the anisotropic model and a constant cou- 
pling coefficient for all the field modes to perform the 
numerical calculation. 

For three dimensional photonic bandgap crystals, if the 
wavevector that corresponds to the bandedge is kg, then 
the dispersion relation in the anisotropic model is de- 
scribed by the following form: — ajc+^|k— kpp, where 
A is a model dependent constant and lOc is the band edge 



frequency. Based on this dispersion relation, the corre- 
sponding DOS is given by: D{uj) = -j^^/cu — a;c0(w — 
ujc)- The memory functions under the anisotropic model 
also can be derived as: 

G(c.) = P'/' Z^^ = (14) 

G^icu) = (3^' ^■^ = (15) 

where /3^/^ = q^~~^372 Vj ^^^d we have used the space av- 
erage coupling strength 77 = ^ J diljci-Epin the deriva- 
tion. 

From Fig. ^ we can see that the spectrum of the 
memory function G{uj) are non-uniform and asymmetric 
as we expect. When the frequency is below the band- 
edge frequency Wc, the memory function becomes pure 
imaginary, indicating the inhibition of the spontaneous 
emission inside the bandgap. The spectrum for another 
memory function Gc{uj) is also similar. It can be eas- 
ily checked that the full-width-half-maximum (FWHM) 
bandwidth of the memory functions in Eq. (|14|l and Eq. 
I|15|l are Aujc- For the bandgap in optical domain, the 
order of ujc is about 10^^~^^ Hz, and the typical lifetime 
of the atom is from 10~^ sec to 10^^ sec, which is much 
longer than the response time of the memory functions. 
Therefore it should be possible to approximate the two- 
time operator products in Eqs. H3I5I) by the equal time 
operator products with the introduction of the Liouville 
operator expansion to be given below. 

For a two-level atom system described by the Hamil- 
tonian H , the time evolution of the atomic operators can 
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be written in general as: 



.(0 = E 



Tl = 



where the Liouville operator C is defined as 



1 



(16) 



In this paper, we will only consider the case in which the 
atom is with a longer lifetime and the pumping is not 
extremely high (i.e., the inverse of the Rabi frequency 
is much larger than the time scale of the memory func- 
tions). In this way we can be sure that the time scale 
of the atomic evolution will be always much longer than 
the time scale of the memory functions. Under such as- 
sumptions, it should be legitimate for us to simply apply 
the zero-th order perturbation terms. This is equivalent 
to use the equal time operator products to replace the 
two-time operator products. We have also checked the 
accuracy of the results by including the first order pertur- 
bation terms. The formulation is given in the Appendix 
and the numerical results only show difference when the 
intensity pumping is extremely high. 

It should be noted that the approximation we have 
made is valid only when the time scale of the system 
response is much slower compared to the time scale of 
the memory function. However, with the use of our ap- 
proximation, we still keep the finite response time of the 
memory function in the formulation. This allows us to 
study the effects that are not considered in the usual 



Markovian regime, where the memory function is simply 
approximated by a delta-function. Although our approx- 
imation includes only one portion of the non-Markovian 
nature of the problem, it should still be quite valid for 
the fluorescence spectrum calculation considered in the 
present work since here the memory function time scale 
is typically much shorter than the atomic response time 
scale. 

Under the zero-th order Liouville operator expansion 
approximation, <7z{t) « azit'), and cr±{t) w cr±(i'), and 
with the identities of Pauli matrices, the generalized op- 
tical Bloch equations in Eqs. (|3I5|I can be reduced to: 



dt'G{t-t')a^{t') +n^{t) 



— oo 

t 



dt'Gc{t-t')a+{t')+n+{t) 

— OO 



dt'[G{t - t') + Gcit - t')]il + a,{t')) + n,{t) 



Please note that these equations are now in a linear form 
with the non-Markovian memory functions. By using 
Fourier transform, we can directly solve these modified 
optical Bloch equations as follows: 



Miuj) ■ X{uj) = A'o(w) 



(17) 



where 



J 



'i{uj + A)+ G{uj) 
Miuj) = ( 

-in 





-i{uj - A) + Gc(w) 

in 



d^{uj + A) 
X{uj)=\ o-+(tj-A) |, and A'o(w) 
(fz{uj) 



i^ 

_ 2 



-iLu + G{uj) + Gc{uj) 

n-iu) + A) 
fi+(ijj — A) 
2n[G{Lu) + Gc{i^)]S{uj) + n^(w) 



r 



and n^{uj), n_|_(ti>), riz{i^), G{w), and G'c(w) are Fourier spectively. The solutions of Eq. II17|I are 
transforms of n^{t), n+{t), nz{t), G{t), and Gc{t), re- 



o-_ (lu + A) 
~ A) 
d-z (w) 



{2g^h^ + »^)n'_^ + n'^n'+^ + i g^nfizjuj) - i2TTg^nG'^5{uj) 
n^n'.^ + (2/^/1^ + n'')n'+^ -if^nfizjuj) + i2^fMG'J{i^) 
2ig^nn'^^ - 2i f^nn'+^ + 2f^g^nz{uj) - ^tt f^g^jG' ujK^) 



(18) 
(19) 
(20) 
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Here we have used the following shorthand notations: In this way the fluorescence spectrum can be easily de- 
termined from Eqs. H18I19|I after determining the noise 

fui = f{^) = — zA + G{uj) correlation functions. It should be noted that here we 

9uj = g{^) = —ioj + i^ + Gc{iLj) cannot directly apply the quantum regression theorem 

^ /if)— ■ G( ) G ( ) since it is invalid for non-Mar kovian process. We avoid 

^ — [uj) — luj + [uj) + c\^) ^Yiis difficulty by linearize the Bloch equations with the 

n'±u — n'±{Lo) = h±(Lo ^ A) Liouville operator expansion and by solving the linearized 

G' — G'{lo) = G{lo) + G iuj) equations directly in the frequency domain. 



IV. FLUORESCENCE SPECTRUM 

Because the two-time correlation function of the 
atomic dipole is proportional to the first order coher- 
ence function ^'^■'(t) [13 of the radiated photon field and 
the fluorescence spectrum can be obtained by taking the 
Fourier transform of the first order coherence function, 
one has: 

/oo 
dTg^^\T)e^-' oc {a+{^ya^{-Lu))R (21) 
-00 

I 



As a check, we first use our formulation to calculate the 
results for the free space case, where the memory func- 
tions are delta functions, i.e. l^fepe^*^** — T5{t) and 
r is the decay rate of the excited atom. The correlation 
functions of the noise operators at zero temperature are 
also delta- function correlated (white noises). Therefore, 
the fluorescence spectrum at steady-state is given by: 



((T+(w- A)(t1(-cj + A))7? = 



A2 



5{uj + A) 



(22) 



r2 + (u; + A)2) 



2A{V^[A - 2{uj + A)2]2 + (cj + A)2[172 A2 |r2 -{UJ + A)2]2} 



where A = ^^ + IS?' + In the limit of strong on- 
resonance pumping (SI ^ F, A = 0), Eq. can be 

reduced to: 

{a~+{u)a-{^^))R = 27r • [^^^5{uj) + (23) 

3 —Y^ 3 

l§± I 4^ I 1 

{oj + r!)2 + Ar2 + 1 r2 ^ + Ar2 ^ 

So far we have proved that the resonance fluorescence 
spectrum exhibits the MoUow's triplet spectral shape 
[l2|: three Lorentzian profiles with peaks in the ratio 
1:3:1, and widths of |F, F, and |F. There is a 
contribution from the elastic Rayleigh scattering in the 



center part (which is a delta function with zero detuning 
frequency) and a contribution from the inelastic Raman 
scattering (the three peaked profiles). The line- width of 
each peak is proportional to the decay-rate of atom and 
the separation of adjacent peaks is proportional to the 
Rabi frequency. This check provides a good support for 
our new formulation. 

Next we come back to calculate the case for a pho- 
tonic bandgap material. By using the density of states 
for three-dimensional photonic bandgap crystals, D{lo) = 
^3/2 V'^ ~ '^cQ('^ — ^c), we get the following noise corre- 
lation functions in the frequency domain at zero temper- 
ature: 



(ri_(cji)ri+(-a;2))i?, = ■nN{u)i)Q{ui + uja - iJc)5{i^i - UJ2) (24) 

{nz{L^i)nz{~^2))R = A^(wi)[27r(5(wi - W2) + (^^(wi - W2))ii]6(wi - Wc) (25) 

(n,(a;i)n_(-a;2))i? = (26) 

{h_{uJl)hz{~UJ2)) R = N{uJi){d^(uJi~ UJ2))RQ{uJl+UJa- ^c) (27) 

{nz{uji)n+{-uj2)) R = N{uji){a+{uji-L02))RQ{'^i+i^L + -i^c) (28) 

{fl+{ui)fl,{-UJ2))R = (29) 



6 




where N{lu) = 4/j3/2 V^a+^- i£e ^ 

After applying these noise correlation functions in Eqs. 
(I24l29f) as well as the memory functions in Eqs. 114t 
I15|l , we can get the fluorescence spectrum under resonant 



FIG. 3: Fluorescence spectra for uja near the band edge ujc- 
Solid line: uja = 100.27/3. Dotted line: uja = 1000/3. The inset 
is the enlarged profile of the lower frequency peak. 



(,/o + 9o) + 2jQ9ohor 



By using Eq. (jSOJ, in Fig. |21we plot the resonance flu- 
orescence spectrum at a constant Rabi frequency when 
the atomic transition frequency tOa is above the band- 
edge u)c- Here the linewidth of each peak is proportional 
to the normalized frequency /3 defined bellow Ea. p5|) 
and the separation of each adjacent peaks is also deter- 
mined by the Rabi frequency. The parameters used in 
the calculation are also labeled in the figure. When the 
atomic transition frequency is far away from the band 
edge {tOa S> ijJc), the normal resonance fluorescence spec- 
trum of MoUow's triplets is obtained just as expected. As 
the atomic transition frequency moving toward the band 
edge, the profiles of incoherence scattering processes be- 
come sharper and sharper as there are fewer and fewer 
DOS available. The profile in the lower frequency is not 
only suppressed but also asymmetrical due to the exis- 
tence of the bandgap as shown in Fig. 01 Its residual 
profile exhibits a sharp edge as shown in the inset of 
the figure. It should also be noted that the peak in the 
higher frequency is enhanced a lot as can be clearly seen 
in the figure. Eventually the peak in the lower frequency 



will be totally suppressed when the atomic transition fre- 
quency is moving more toward the band edge. At this 
time the resonance fiuorescence spectrum now only has 
two peaks, as shown in Fig. 2| This is of course again 
due to the bandgap effect. It is interesting to see that 
now the enhancement of the original middle frequency 
peak becomes smaller than the original higher frequency 
peak. If the atomic transition frequency moves below 
the band edge, then it is no longer possible to resonantly 
pump the two level atom located within the photonic 
bandgap material. Although in principle non-resonant 
pumping through the use of three-level atoms can be real- 
ized, the fluorescence spectra will be different when com- 
pared to the resonance fluorescence spectra considered 
in the present paper. This will be one of the interesting 
topics that can be further studied in the future along this 
line of research directions. 

We also plot the total resonance fluorescence power 
spectrum in Fig. |31 by integrating over all frequencies. 
The curve exhibits three regions of different behaviors for 
different offsets between the transition frequency and the 
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FIG. 4: Resonance fluorescence spectra for uja smaller than 
the bandedge u>c at a constant Rabi frequency. 
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FIG. 5: Total fluorescence power spectrum for cua near the 
bandedge. 



bandedge frequency. When the transition frequency is far 
above the band edge, Region 1, the total power spectrum 
remains almost the same due to small variation of avail- 
able DOS. As one moves the transition frequency closer 
to the band edge, Region 2, the total power spectrum 
will increase as the photonic crystals begin to impact the 
inelastic scattering processes, such as the case in Fig. O 
The total power spectrum will saturate when the first 
profile of MoUow's triplet is totally suppressed, Region 
3. For comparison, the optical power corresponding to 
each profile peak is also plotted in Fig. for different 
wavelength offsets. 



CONCLUSIONS 
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In conclusion, we have developed a new formulation to 
calculate the fluorescence spectrum with non-Markovian 
photon-atom interactions, and have successfully applied 
it to the case of a single atom embedded in a photonic 
bandgap crystal. By introducing Liouville operator ex- 
pansion, we can overcome some of the difficulties caused 
by the non-Markovian nature of the problem due to the 
non-uniform distribution of the photon states in the pho- 
tonic bandgap crystals. Although the approximation we 
have made is valid only when the time scale of the sys- 
tem response is much slower compared to the time scale 
of the memory function, however, with the use of our ap- 
proximation we can still keep the finite response time of 
the memory function in the formulation. This allows us 
to study the effects that arc not considered in the usual 
Markovian regime, where the memory function is simply 
approximated by a delta- function. Although our approx- 
imation includes only one portion of the non-Markovian 
nature of the problem, it should still be quite valid for 
the fluorescence spectrum calculation considered in the 
present work since here the memory function time scale 



FIG. 6: Individual fluorescence powers for different peaks. 
Solid line: power of the lower frequency peak; Dashed line: 
power of the central frequency peak; Dotted line: power of 
the higher frequency peak. 



is typically much shorter than the atomic response time 
scale. We think one important contribution of the present 
work is to demonstrate that even in this unconsidered 
not-fuUy-non-Markovian regime there are still some in- 
teresting new phenomena like the suppression and en- 
hancement of the relative fluorescence peak amplitudes 
at different wavelength offsets. Historically, in 1975 the 
MoUow's triplets for the free space case was observed 
experimentally for the first time by using the hyperfine 
structure of the sodium atom It shall be very in- 

teresting to see if one can actually observe the predicted 
phenomena experimentally in the near future. We also 
believe that the results we have found should open up 
many new studies along this line of research directions. 
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APPENDIX: FIRST ORDER LIOUVILLE 
EXPANSION 

In this appendix we carry out the Liouville operator 
expansion in Eq. H16() to the first order: 

H^Ho = ^((T_e^^* + a+e-'^*) + ft^c^fe4afe (A.l) 

k 

I 

+ 



It can easily be seen that above equations are reduced 
to the zero-th order equations when the Rabi frequency 
is small. With the first order Liouville operator ex- 

I 



The atom operators (Tz(t), cr-(t), and cr+(i) can then be 
approximated by 



(A.2) 
(A.3) 

(A.4) 

I 

pansion, the generalized optical Bloch equations in Eqs. 
(|3I5II become: 



coa[n{t - t')]a,{t') + ism[n{t - i')][CT-(t')e'^*' " t^+(i')e"'^*'] 

cos[ri(t - t')] [a_ (t')e'^' - {t')e-'^' 1 + 2 (^') + ^"^'^^ (^')] 
sin[f^(i-t')]<T,(t') 
^cos[17(t-t')][a_(t')e'^' -a+(t')e-^^' ] + - + e^'^* 



CT-{t)=i-a,e-'^' - f dt'G{t-t')cos[n(t-t')]cr^(t')-i f dt'G{t - t') sinn{t - t')]{^-^-^^) (A.5) 

cr+(i) = -i-cr,e*^* - f dt'Gc{t~t')cos[n{t-t')]a+{t')+i f dt'Gc{t - t') smn{t - t')]{^^^-^^){A.6) 
2 J-00 J -00 2 

d,{t)=in{G^e'^' -a+e-'^')- f dt'[G{t - i') + Gc(t - t')]{ ^ ^ ^"^^'^ ) (A.7) 
/" dt'sin[f^(t-t'){G(i-^')'^-(i')e''^*' -Gc(t-t')'T+(i')e"''^*'} 

^ — 00 

- f dt' cos[r!(t - t')][G{t - t') + Gc(t - i')](l±^il^) 



By the same method of Fourier transform, we can get lation functions of the noise operators at zero tempera- 
the fluorescence spectrum with following non-zero corre- ture: 
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(n_(cJi)n+(-cj2))i?, 

{h-{uJi)nz{-uj2))R 
(n^(wi)n+(-W2))fl 



= Ni{uji)[2Ti5{uji - UJ2) + + (f?'+(^i-c.2)>fl]0('^i +i^a + ^-^c) (A.8) 

+ A^2(wi)[27r(5(wi - tJ2) - {(^' -(UJ1-UJ2)) R ~ (f^' + (c^i-t^2))fl]0('^l +c^a - - t^c) 

= Nl{uJl)[2Tl5{uJl ~ LO2) + 2(CT'_(^,_^,))fl + {a,{uJl - UJ2))r]Q{uJi +UJa + ^-LOc) (A.9) 

+ N2{LOi)[2'n:5{LOi - tj2) - "^{c^' -(^1-^02)) R + (^^2(^1 - t^2))i?]B(wi + - O - tJc) 

UJa + CJi 

= Ni{uJi)[2Tl5{uJi - tJ2 - A) + 2(CT_(a)i - UJ2))r + (^7(^,_^,))i?]©K + t^a + - Wc)(A.10) 

+ N2{ijJi)[-2ti5{uJi - tJ2 + A) + 2(ct_(cJi - UJ2))r - (^7(tui-i^2))-R]Q(^l + t^a - - Wc) 

= Ni{uj^)[2t:5{uj^ - t^2 + A) + 2(a+(a)i - tj2)>i? + (^;(^,-^,))fl]0('^2 + + - ^c)(A.ll) 

+ N2{uji)[-2ti5{uji - W2 + A) + 2(ct+(wi - W2))/? - (^|(^i-c^2)>fl]0('^2 + cja - - t^c) 



where 

(jJa + W + ii 

cr'±(wi-t^2) = o-± (tJi - tj2 T A) 

^From Eq. ljA.8IA.lip , the correction of high order terms 
wiU have impacts on the fluorescence spectrum only when 



the Rabi frequency f2 is of the same magnitude with other 
characteristic frequencies of the system, i.e. LOa, and lOc- 
For optical systems we are interested in, the zero-th or- 
der Liouville operator expansion shall give us reasonable 
results. Only for the cases of very strong pumping power 
or ultra short atom lifetime we need to take into account 
the higher order corrections. 
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